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Abstract 

We study the problem of triviahty in the four dimensional Nambu-Jona-Lasinio model 
with discrete chiral symmetry using both large- expansions and lattice simulations. We 
find that logarithmic corrections to scaling appear in the equation of state as predicted by 
the large-A^ expansion. The data from 16^ lattice simulations is sufficiently accurate to 
distinguish logarithmically trivial scaling from power law scaling. Simulations on different 
lattice sizes reveal an interesting interplay of finite size effects and triviality. We argue that 
such effects are qualitatively different for theories based on fundamental scalar rather than 
fermion fields. Several lessons learned here can be applied to simulations and analyses of 
more challenging field theories. 
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1. Introduction 



The Nambu- Jona-Lasinio (NJL) model [1] describes the dynamics of massless fermions 
that have a contact interaction. The theory has two phases. At weak couphngs fermions 
are massless, while at strong couplings they acquire dynamical mass due to spontaneous 
chiral symmetry breaking. The transition from a massless to a massive phase is accom- 
panied by the appearance of light scalar (and pseudoscalar) bound states. Through the 
exchange of these composites an effective interaction develops over physical scales and the 
cutoff disappears from the model [1]. The original theory, although not renormalizable 
in perturbation theory, becomes renormalizable and appears to have the same low energy 
limit as a linear u-model [2]. Since the effective interaction is of the Yukawa type, it is 
important that, measured in physical units, the mesons retain a finite extent in order to 
ensure an interacting low energy limit. This is possible to achieve below four dimension 
[3,4], but for d > 4, the mesons are always pointlike in the continuum limit and the cor- 
responding theory is free. In the marginal dimension, (i = 4, there are calculations and 
analyses supporting triviality as well [3,5], but a rigorous demonstration is lacking. 

The reason for triviality in the NJL model is different from that in the more familiar 
and more extensively studied case of 0^ theory [6] . The essence of the chiral phase transi- 
tion is to tune the coupling in such a way to allow for the formation of bound states. Once 
bound states appear in the spectrum, everything else related to chiral symmetry breaking 
follows naturally [7]. Keeping the size of the bound state finite, on the other hand, is 
a matter of balancing the attraction due to the interaction and the zero-point repulsion 
due to the kinetic energy. Above four dimension, this is not possible to achieve since the 
short-distance attraction becomes progressively stronger with increasing dimensionality 
and collapsing bound states always lower the energy. 

In four dimension, however, the attraction and the repulsion scale identically and 
the physics is driven by whatever is left over after their cancellation. So, the problem of 
triviality in four dimension amounts to more than simple power counting and a decisive 
answer requires a careful study. The technical difficulties lie in the fact that power counting 
in the four dimensional case is modified by logarithmic corrections and these logarithms 
are notoriously difficult to establish both theoretically and numerically. There are several 
techniques, such as large- AT, e-expansions or perturbation theory, that can provide some 
hints of how the logarithms enter the game [3,8,9]. In the NJL model the physics of 
triviality is intrinsically nonperturbative due to the presence of bound states and chiral 
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symmetry breaking - the weak coupling phase is symmetric and the Goldstone physics is 
not accessible to perturbation theory. Thus, the only analytic method that can be of some 
utility in the context of NJL models is the 1 /N expansion which we will use as a guide in 
analyzing the data. 

In the past the NJL model has captured considerable attention and has stimulated 
activities in a wide range of contexts, from purely theoretical [3,10-12] to practical model 
building, such as nuclear physics and the intermediate energy sector of QCD [13,14], 
technicolor [15] and the top condensate model [16], to mention a few. Most of the studies 
have been performed in the large- A?" limit, for the reasons mentioned above, and systematic 
studies of the 1/N corrections have been carried out. The 1/N expansion indicates that 
below four dimension the theory is renormalizable [3,8], whereas for d > 4 it is trivial. 
There is little doubt that this changes at low- A?", but no proof outside of the 1 /N expansion 
exists. In that sense we have undertaken this project in order to gain some insight and 
develop some criteria to find logarithms in a theory that is known to have them. 

Recently several groups have carried out numerical simulations of various versions of 
the NJL model, at [16] and below four dimension [18,19], but the problem of triviality in 
four dimension has not been tackled intensively before. In this paper we want to address 
precisely this issue with an emphasis on how the logarithms should be tracked down. We 
should point out that we do not expect any surprises as far as the large-A/" predictions 
are concerned and, in that sense, we anticipate the outcome that the theory is trivial. 
Our main motive is to see, knowing the answer, with what degree of confidence we can 
establish it numerically. Our main efforts will be concentrated on the equation of state 
(EOS) and the logarithms that accompany critical scaling. We have been encouraged 
to undertake this project by the recent observation that the logarithms in models with 
composite mesons, (therefore the NJL model), show up in a qualitatively different way than 
in theories with elementary scalars, like [20]. Interestingly, the logarithmic behavior 
from our simulations on large lattices (16"^) agrees with the composite meson case, as 
expected. However the logarithmic behavior from the simulations on the smaller lattices 
(8^) is consistent with the elementary scalar case. This apparent disagreement between 
two lattice simulation results must be due to the systematics of finite size effects and 
understanding the reason behind it will increase our confidence in deciding whether a 
particular lattice result is tainted by finite size effects. Our efforts include some technical 
aspects of this problem as well. In particular, we will show how to distinguish power 
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law scaling from mean field scaling modified by logarithmic corrections. This problem is 
"reversed" from searches for a nontrivial fixed point in other models. As will be discussed 
in section 4, when we discuss the lattice data, the logarithms in the NJL model reveal 
themselves in a clear and unambiguous way. Especially interesting, in this context, is the 
manner in which the power law hypothesis fails. 

There are important lessons to be learned from this exercise alone. They concern 
numerical simulations of theories with nontrivial fixed points for which no analytic tools 
arc available. On a technical level, establishing the logarithms (therefore triviality), is a 
more demanding job than diagnosing their absence in theories that are not trivial. 

Setting our goal as "logarithm hunting" , we have chosen to simplify our task as much 
as possible. For that reason, we decided to work with N = 12 flavors in order to remain 
close to the large- A?" limit. This turned out to be a profitable decision for another reason. 
Tunneling is suppressed at large-A, so we could work relatively near the model's critical 
point without suffering from this particular finite size effect. Also, the choice of discrete 
chiral symmetry, instead of continuous chiral symmetry enables us to work in the chiral 
limit and avoid the extrapolations to zero bare mass which complicate other studies in the 
field. 

The layout of this paper is as follows. In Sec. 2 and 3 we solve the model in leading 
large- A" order and identify the logarithms of triviality in the equation of state, suscepti- 
bility, and specific heat. These results are contrasted to the analogous quantities in cp'^. In 
Sec. 4 we present the simulation data and fits on 16"^ and 8^ lattices. Wc shall see that the 
16^ lattice simulations expose the logarithms of the largc-A expansion with considerable 
precision. The 8^ data is strongly distorted by finite size effects which in some cases ob- 
scures the logarithms in curious ways. In Sec. 5 we discuss our results and the lessons we 
have learned which can be applied to studies of more challenging field theories. 

2. Large-A limit 

In this and latter sections we review the large- A limit of the NJL model. The emphasis 
will be on the scaling violations and how they enter the EOS and expressions for other 
thermodynamic quantities. In order to facilitate the large-A expansion, the four-fermi 
lagrangian can be written in the form of a Yukawa theory 

L = il;{-i^ + m + ga)ilj-^a'^ (2.1) 
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The original, four-fermi, lagrangian is recovered after integrating over the a field. Since 
L is quadratic in fermionic variables, we can integrate them out and obtain the effective 
action 



S^j^ = 1 / _ iVtr ln(-z^ + m + £/(t) (2.2) 

In the large- A'^ limit the dominant configurations are obtained by the saddle point method. 
The resulting equation of state (EOS) is the minimum condition a — Ngtr J^l/{^+m + 
ga) = 0. After identifying the fermion mass as S = m + ga, the minimum condition reads 

E = m-g^ <i;ip> (2.3) 
The curvature of the effective action is the inverse susceptibility 
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The appearance of soft modes at the critical point is manifested in the divergence of the 
susceptibility and the vanishing of the fermion mass. Setting E = in eq.(2.4) gives 
1 — Ag'^ J = 0. With this expression for gc, the EOS can be recast in the form 



where t = g'^/g'^ — 1. This is a convenient form for extracting the critical indices since one 
need only count the infrared divergences of the right hand side. The indices associated 
with the EOS are defined as [21]: < ipip >~ and x ~ t~'^ at m = 0, and < ipip >~ m^^^ 
at t = 0. Since there are only two fields the system responds to, only two out of three 
indices are independent. They are related by the scaling relation P{5 — 1) = 7. Above 
four dimension, the integral in (2.5) is infrared finite, and the exponents do not change 
with dimensionality - they have the mean field values: /3 = 1/2,(5 = 3,7=1. Below four 
dimension the right hand side scales as E''"^ resulting in nongaussian critical behavior with 
indices /3=l/(d — 2),5=ci— 1,7 = 1. In four dimension, the integral is logarithmically 
divergent giving rise to scaling violations. The EOS in this case reads 

-+t = a/s2ln— (2.6) 

where a and b are constants, and the cutoff is set equal to one. Using a momentum cutoff 
gives a = l/47r^ and 6=1. In the chiral limit, the fermion mass scales as the order 
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parameter, S ~< ipip >. At the critical point E = (m + Cm}/^), but since 5 > 1, the 
linear term is subleading and the scaling of S and < ipip > coincide again. The two limits 
of the gap equation give 

t = ag^'S^ In ^ (m = 0) (2.7a) 

Zj 

m = ag^'S^ In (t = 0) (2.76) 

Zj 

Up to logarithms the exponents have their gaussian values (3=1/2 and 5 = 3. 

There are two universal quantities to be extracted from the susceptibility. In the 
critical region, in the chiral limit, the susceptibility scales as x ~ C\t\~'^ . The amplitude 
C is different in the two phases and the ratio C+/C_ is a universal quantity. Mean field 
theory predicts C_/C+ = 2 [21]. In the symmetric phase it is easy to compute both C_ 
and 7. Setting S = in eq.(2.4) and using the definition of the critical coupling, we get 



X 



-' = l-9V9l-{-t) (2.8a) 



which implies 7 = 1, C_ = 1. In the broken phase the first two terms in eq.(2.4) cancel 
because of the gap equation. Combining the third term in eq.(2.4) with the gap equation 
(2.5), gives — 2t ~ S^and the expression for the susceptibility becomes 



Clearly, the exponent 7 = 1 is the same as in the symmetric phase. There are no log- 
arithmic corrections to scaling for the susceptibility. The amplitude, however, receives 
corrections due to scaling violations so that 

Cl^'^~ ln(6/S2) (^-^^ 

The wavefunction renormalization constant is given by the k'^ term of the a propa- 
gator. Due to radiative corrections such a term is generated in the 1/N expansion. The 
leading log contribution is given by 

The anomalous dimension r/, defined by Z ~ S'', vanishes up to logarithms. 
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There are basically two candidates for the correlation length in this model: the fermion 
and cr- masses. It is simple to verify that they both scale the same way. The mass of the 
a particle, for example, is obtained through the second moment %• Using eqs.(2.8) and 
(2.10) leads to = Zx~^ ~ The scaling of the fermion mass, S, is obtained easily 
from the gap equation, (2.7a). It follows that i' = and so v = 1/2 with the same scaling 
violations as for the order parameter. 

The specific heat is obtained from the effective action, eq.(2.2). Like the susceptibility, 
it too is characterized by two quantities, the exponent a and amplitude A. In the critical 
region its scaling is given by C = A\t\~°'. The amplitudes are different in the two phases 
and their ratio A-/A+ is universal [21]. Mean field theory predicts o; = 0, A- — and 
the specific heat undergoes a jump. Unlike other exponents, a can be either positive or 
negative. In four dimension logarithmic corrections affect the scaling of C. However, to 
leading order in A_ remains zero. It is easy to see this since in the symmetric phase, 
< a >= 0, and the action is independent of g. Thus, A_ = 0. In the broken phase, 
straightforward algebra leads to 

C — (2.11) 

ln(l/t) ^ ^ 

Higher order corrections might change the power of the logarithm, but it remains in the 
denominator. For scalar theories there are several results on the amplitude ratio and 
critical index a (e-expansion and large- A^) [22]. In the symmetric phase A_ = 0{1/N), 
so it does not appear in leading order. The amplitude ratio, therefore, scales as 1/N. 
In four dimension, ct = as well, but with logarithmic corrections of the form C ~ 
l/{lnt)^^~^'^^^^'^^\ Thus, depending on A^, the power of the logarithm can change it's 
sign and the specific heat can change its behavior from vanishing to divergent. The e- 
expansion studies indicate that this does not happen in Yukawa models [19]. 

3. Triviality of the four-dimensional theory in the large- limit 

To establish the connection between scaling violations and triviality, we introduce the 
renormalized coupling. It is a dimensionless low-energy quantity that contains information 
about the non-gaussian character of the theory. It is conventionally defined as [23] 

(nl) 
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where the nonhnear susceptibihty x''"'^^ is the zero-momentum projection of the connected 
four-point function 



X("') = ^ = < cT(0)a(l)a(2)a(3) (3.2) 

The normahzation factors, X = < cr{0)a{x) >c and in eq.(3.1) take care of the four 
fields and the three integrations. In a gaussian theory all higher-point functions factorize, 
so gn vanishes. Using the hyperscaling hypothesis, this can be converted into. 



9R ~ e^2A-7-d.)/. (3 3) 

where A = + 7 is the gap exponent. Being dimensionless, gn should be independent of 
^ if ^ is the only scale. Thus, the validity of hyperscaling requires that the exponent must 
vanish. It implies the relation, 2A — 7 — di^ = 0, between the critical indices. In general, 
it is known that the following inequality [24] holds 



2A < 7 + dz/ (3.4) 

The exponent in the expression for gji is always non-positive, so that violations of hy- 
perscaling imply that the resulting theory is non-interacting. Above four dimension, the 
exponents are gaussian (7 = 1, A = 3/2, = 1/2). In this case, it is easy to verify the 
above inequality: 3 < 1 + d/2, which amounts to d > 4. In four dimension there are 
logarithmic corrections to scaling that are believed to drive gn to zero. Scaling violations 
in any thermodynamic quantity propagate into the renormalized coupling and, according 
to eq.(3.4), these violations lead to triviality. Because this is a marginal case, it requires 
special treatment. 

In order to make our point simple, we illustrate how logarithms appear in the effective 
actions of two solvable models: 0^ and {'ipip)^ theories both in the large-A?" limit. The 
effective actions for the two models are [5,25] 

1 m"^ a 

V{a) = —ta^ + a^ln{l/a) (3.56) 

This is the leading log contribution only. In the first example, it is clear how log-corrections 
lead to triviality. The logarithm can simply be thought of as coming from the running 
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coupling - quantum corrections lead to the replacement A — > Ajj. The vanishing of the 
renormalized coupling is then manifest from eq.(3.5a). In the case of fermions, eq.(3.5b), 
the details are completely different - the analogous reasoning would lead to an erroneous 
conclusion that the renormalized coupling increases in the infrared. This form is obtained 
from the effective action, 5'e//(cr), in eq.(2.3). It is easy to see how the logarithm in (3.5b) 
appears in the numerator. For example, after taking the second derivative of x~^, eq. (2.4), 
the last term generates a logarithmically divergent contribution which is the log-term in 
eq.(3.5b). The vanishing of the renormalized coupling here follows from the wave function 
renormalization constant Z ~ l/ln(l/(T) [5]. 

In both cases the renormalized coupling is obtained through the nonlinear suscepti- 
bility. The correlation length is related to the susceptibility by = x/Z. For magnets 
(i.e. scalar theories like 0^) the following relations hold, 

~ ^'riTm) ~ WW) '■^■'^''^ 

so gn vanishes at a logarithmic rate. 

For fermions, on the other hand, we have 

X^-'^ ^ ln(l/^), Z ~ — ^ (3.7a) 

ln(l/cr) 

gj,^ZHn{l/a)r.:^-lj-^ (3.76) 

and again gj^ vanishes at a logarithmic rate. In this context the following point should be 
made. The nonlinear susceptibility in the NJL model is a connected four-point function 
for the composite ijjijj field. The free fermionic theory is not gaussian in so even in 
free field theory gn, defined in this way, does not vanish. The fact that gR ^ near the 
critical point indicates not only that the resulting theory is gaussian, but that it is also 
a purely bosonic theory. So, for the NJL model, the bare theory and its continuum limit 
bear no resemblance. This is quite different from what happens in the (p"^ model where 
the resulting continuum limit has the same physical content as the bare theory but with 
vanishing coupling constant. 

The important point in this comparison of the two models is the fact that the loga- 
rithms in the EOS appear in different places. This property is generic for the two models 
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and persists beyond the 1/N expansion. In order to emphasize this point, we consider the 
critical EOS for both models. They are obtained from the effective potential by simple 
differentiation. To make the connection with 5, we take t = 0. The critical EOS for the 
magnets is [25] 

This defines the exponent 5. Because of the scaling violations, eq.(3.8) vanishes faster 
than a pure power. So the "effective" 5 is bigger then its mean field value. It is easy to 
understand the correction to scaling in eq.(3.8). 1/5 is simply the logarthmic derivative 
of the order parameter 1/5 = dlnM/dlnh. Since the magnetization saturates at strong 
fields, its slope vanishes in the large-/i limit. Thus, 1/5 decreases away from the h = 
limit. As was shown in ref.[20], such corrections to scaling can never occur in the case of 
the chiral transition. 

For the NJL model the critical EOS reads, 

m - cT^log(l/cT) (3.9) 

Unlike scalar theories, the log's appear in the numerator - the right hand side in eq.(3.9) 
vanishes slower than the pure power and the "effective" 5 is smaller than the (pure) mean 
field value. This is also easy to understand since at large mass < ipip >~ m. Thus, away 
from the chiral limit 1/5 is driven to 1. As pointed out in ref.[20] this difference in the 
position of the logarithm, eqs. (3.8-9) or, equivalently, the sign of the scaling violations 
is generic for the two models and is intimately tied to the fact that in NJL scalars are 
composite and in (f)^ they are elementary. 

The distinction between the two models proved to be useful in analyzing the lattice 
data. In particular, if the NJL model is simulated on a small lattice, then, for a fixed 
mass, the right hand side of eq.(3.9) will be smaller than it is in the thermodynamic limit 
because finite volume effects tend to restore the symmetry, and reduce the order par am- 
eter. If the infinite volume a is to be fitted to the data in the form m ^ a , then the 
appropriate exponent 5', which might fit the data, must be bigger than its actual value. 
This larger value of 5 could also be mimicked by mean field scaling supplemented with 
the log in the denominator. Such data should be not be trusted since their behavior is 
qualitatively different from the behavior of the system in the thermodynamic limit. As 
the lattice size is increased, the effective 5 would eventually fall below its mean field value 
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suggesting the proximity of the thermodynamic hmit. Such data, although not necessarily 
completely free of finite size distortions, do share the same qualtiative features with the 
physics of the infinite volume. So, for fermions the bound 5 < 3 separates the data into 
two groups and tells us which one is more reliable. This is possible because finite size 
distortions and the approach to the true scaling law have the opposite tendencies. For 
scalar theories, on the other hand, the two have the same tendency. Nontrivial scalar 
theories have 5 > 3, to begin with. For example, from eq.(3.8) it is apparent that in case 
of logarithmic corrections to mean field scaling, the finite size distortions which tend to 
reduce the magnetization can be masked by adjusting either the exponent 6 or the strength 
of the logarithm so that the right hand side of the critical EOS assumes an appropriate 
magnitude compatible with the given h. In fact, strong finite size effects would force the 
fit with 5 > 3 even if the scaling is mean field, so work on a larger lattice is necessary 
to either confirm or disprove this result. Thus, in scalar theories a careful finite size scal- 
ing analysis seems to be necessary to disentangle the log's from the finite size effects. In 
that sense, establishing triviality in 0"^ is more delicate than it is in fermionic theories [6] . 

4. Simulation results 

We carried out extensive simulations of the model described by Eq.(l) on 8^ and 
16^ lattices. The small lattice simulations were done to monitor finite size effects. Of 
course, many more calculations on a variety of lattice sizes could be done to make our 
considerations on finite size effects more quantitative, but 8^ and 16^ simulations will 
prove adequate for our purposes. If Four Fermi models continue to play an ever increasing 
role in high energy theory, then more thorough studies would be needed. 

Although there are some cases in which the finite size effects were large and curious as 
we shall see below, most of the 16^ simulations were relatively free from finite size problems 
for rather wide ranges of couplings and yielded precise continuum scaling laws. In fact, the 
logarithms of triviality will be seen clearly in our study of the Equation of State. There are 
several reasons for this success. First, the model can be simulated for zero quark mass. So 
we are free of the difficult issues of extrapolations to the chiral limit. Second, the discrete 
character of the symmetry breaking made the Hybrid Monte Carlo algorithm particularly 
effective and allowed us to run the simulations with a relatively large time step (dt=0.1) 
while still having a high acceptance rate (from 60 to 90 percent was typical). Third, we 
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could dial the number of flavors large to compare to the 1/N expansion (we chose Nf = 12, 
as in our past work in three dimension [18]) which suppresses fluctuations enormously. 

As discussed in the theory sections above, scalar systems and fermion systems are 
clearly distinguished by the positions of the scale breaking logarithms in their critical 
behaviors. One of the more interesting features we discovered here is that those distinct 
behaviors can be blurred by finite size effects. Comparison between 16^ result and 8^ 
clearly illustrates such effects. 

The lattice action, the Hybrid Monte Carlo algorithm and other technical details of 
the simulations follow ref[18] in detail. The interested computational physicist should 
consult that reference for details. In our discussion and Tables, we shall add any technical 
information which was special to working in four dimension. 

For clarity, we shall organize our lattice simulation results into various subtopics, 
following the theoretical discussion above. 

Equation of State 

In Table I we have collected the 16^ data for the average action, its specific heat, 
the vacuum expectation value of the a field, and its susceptibility x- We also tabulate 
there the statistics of each measurement and the statistical errors on each quantity. The 
errors were calculated with considerable care using the usual binning techniques, because 
we are interested in high quality fits to distinguish different functional forms in the model's 
critical behavior. Note that in the immediate vicinity of the critical point, we accumulated 
particularly high statistics (1/3 of a million sweeps) to combat critical slowing down. 

We first considered the order parameter's dependence on coupling in the broken sym- 
metry phase. The logarithmic fit, 

P-Pc = aa^i\n(l/a) + bf (4.1) 

was tried for the twelve data points in the coupling range f3 = .54— .595. A very good fit was 
found (confidence level=70 percent) and the parameters were determined: (3c = .6044(1), 
a = 2.25(2), b = 1.16(5), and p = 1.041(12). The data and fit are shown in Fig. la. 

Perhaps the most interesting feature of this result is the fact that the power of the 
logarithm in Eq.(4.1) is determined with such precision and confidence. The fitting rou- 
tine (standard least square method) searched for the best value of the power and found 
1.041(12). Note that the fit does not just accommodate the immediate vicinity of the 
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critical point : the parameter 'b' in Eq.(4.1) which sets the scale of the logarithm allows 
us to fit both the strong coupling region (where a is large and mean field scaling should 
apply) and the critical region (where fluctuations causes the logarithmic correction to the 
mean field behavior of the model). 

It is interesting to plot the data in a slightly different way to show the importance of 
the logarithm more visually. Consider Fig. 2a where we plot the data as {13 — I3c)/<j'^ vs. 
ln(l/a") for the choice /3c = .6044. The numerical significance of the logarithm and the fact 
that its power is very close to 1.00 are apparent. 

We also attempted equation of state fits which are free of logarithms. For example, if 
the theory had pure power law singularities, then the form. 



should apply in the broken phase. When a least square fit of this form is attempted, the 
iterative procedure ( the search for the best confidence level ) does not converge, although 
high confidence levels are found for fixed fitting parameters. The fitting program finds that 
it can always increase the confidence level by increasing 'a' positively, while increasing 'b' 
negatively and driving l/Pmag to 2 from above. Such a tendency is reminiscent of the 
e-expansion where the limit of d = 4 is marked by the fact that power law scaling is 
converted into logarithms through the formula Inx = lime_»o(a;^ — l)/e. The behavior of 
the fit is compatiable with the weaker logarithmic dependence on a found above. Our 
data is sufficiently accurate and it extends over a wide enough range of P to favor the 
logarithmic equation of state over the pure power law formula. 

In Table II we collect our 8^ measurements which parallel the 16^ measurements 
presented in Table I. First consider the equation of state in the broken phase. Fits of the 
form Eq (4.1) and (4.2) which should have worked in both the strong coupling and critical 
regions were unsuccessful ( the confidence levels were infinitesimal ). However, if we took 
a subset of the data in the critical region alone reasonable fits could, of course, be made. 
For example, a pure power law fit such as. 



could be made over the relatively narrow region from /3 = .57 to .595 with a fair confidence 
level (20 percent) and results: /3c = .607, a = 1.56(9), p = .54(3), and b = .002(13). The 



f3c-f3 = acT^/^-"" + 6(7^ 



(4.2) 



a = a{(3c -py-b 



(4.3) 
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fit is shown in Fig. lb. Note that the fit is consistent with mean field theory, unadorned 
with logarithms. 

To gain more insight into the problems with 8^ simulation data, let us plot it as we 
did the 16^ data in Fig. 2a. The result is shown in Fig. 2b. Comparing with Fig. 2a, it 
appears that the 8^ data suffers from finite size effects as {3 approaches the critical point 
and there are only 'hints' of the true logarithmic corrections in Fig. 2b. This plot shows 
clearly that it would not be sensible to 'blindly' fit 8"^ data to any continuum model scaling 
law, even though the 16^^ data shows continuum behavior over the same range of couplings. 

Susceptibility 

The measurements of the order parameter a allow us to compute its variance and 
construct the model's longitudinal susceptibility The data is given in Table I. First 
consider the broken symmetry side of the transition and a determination of the critical 
index 7 assuming a pure power law singularity. For (5 ranging from .54 through .595, we 
fit the susceptibility data with 

= a(/?e -PV + c (4.4) 

We are particularly interested in fits which produce an estimate of the parameter 'c' con- 
sistent with zero, indicating that fic is a critical point. The least square routine gives: 
a = 9.14(1.12), 7 = 0.967(45), c = 0.002(10) for the choice (3c = .6028. If trial values of 
Pc are chosen larger than .6045 or smaller than .6022, then the parameter 'c' is pushed 
significantly away from zero. However, for all such choices, 7 is found to be .967(45). So 
this result is very robust given a power law hypothesis. The data and the fit are shown in 
Fig. 3a. the confidence level of the fit is 63 percent (chi-squared = 7.13 with nine degrees 
of freedom). 

As discussed in the theoretical sections above, there are weak logarithmic corrections 
to power law scaling of the susceptibility in the broken phase. So, it is interesting to try 
fits of the form, 

X"^=af2-— — 4 V/?c-/5)^ + c (4.5) 

^ V ln(l/c7) + 1.167 

where the scale in the logarithm (i.e. 1.16) has been taken from the fit Eq.(4.1). Again, very 
good fits are found. Requiring that 'c' be compatiable with zero, produces the prediction 
(3c = .6040(8), 7 = 1.054(47), and a = 7.41(95). Although (3c is not predicted as accurately 
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as in the Equation of State fits, the critical index 7 is again found for all choices of /?c to be 
1.054(47). The confidence levels of these fits are typically 55 percent. Clearly the logarithm 
in Eq.(4.5) is not very significant numerically due to the presence of the substantial additive 
term 2. If we had a wider range of a values, the logarithm would play a greater role. The 
logarithms of triviality are much more easily exposed in other observables of this model. 

Now consider the susceptibility on the symmetric side of the transition. No logarithms 
are expected in the scaling law here, so we try a fit of the form Eq.(4.4). We find excellent 
fits (confidence levels of 99 percent) with : o = 5.62(1.65), Pc = .6024(7), and 7 = 1.0(1). 
Again, 7 is compatiable with unity, as expected theoretically at large N. The fit is shown 
in Fig. 3a. 

For the 8^ lattice, we accummulated susceptibility data in both the strong and weak 
coupling phases, and we attempted fits of the form Eq.(4.5) that were successful on the 16^ 
lattice. To avoid vacuum tunnelling we did not study couplings too close to the critical 
point. The ranges of /? = .57- .595 and .62-. 655 satisfied this criterion. The data and 
power law fits are shown in Fig. 3b. In the broken phase the fit gave a critical index 7 
= .97(12), compatiable with theoretical expectations. However, the confidence level of the 
fit is barely acceptible (9.6 percent corresponding to a chi-squared of 6.4 with 3 degrees of 
freedom), and the other parameters in the fit were a = 11.3(3.7) and b = —.029(29) for 
a critical coupling of /3c = .6066. Note that the crude strong coupling data shown in the 
figure at /3 = .54, .55, and .56 were not used in the fit. The weak coupling data and its 
fit are also shown in the figure. The confidence level was very good (93.9 percent) and the 
fit predicted the critical index 7 = 1.21(10) with a = 9.8(2.8) and b = .025(8) for /S^ = 
.6066. So the fitted 7 is two standard deviations from the theoretical prediction. Since the 
16^ fits work beautifully with 7 = 1.05(5), this discrepancy is most probably a finite size 
effect. 

In summary, the 8^ susceptibility study is not quantitatively reliable, in contrast to 
the 16^ results which are very close to the continuum model predictions. 

Equation of State at Criticality 

Another feature of the theory which displays the logarithm of triviality in an accessible 
expression is the equation of state at criticality. In Sec. 3 above we contrasted the behavior 
of scalar vs. fermion models in Eq.(3.8) and Eq. (3.9). To distinguish the two behaviors 
we must simulate the model at f3c in the presence of explicit symmetry breaking. We 
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anticipate that this work will be subject to large finite size effects and will not be as 
quantitative as our other measurements which avoided all the pitfalls of the critical point. 
We shall see that our results will favor Eq.(3.9) over Eq.(3.8), and a comparison with 8^ 
measurements will expose some particularly interesting finite size effects. 

We ran simulations at nonzero fermion mass m at /3 = .6034. The results are compiled 
in Table III. It would be interesting to repeat this work at other estimates of Pc, but 
high statistics are needed since we are working very close to criticality and statistical 
fluctuations are large. We begin by plotting the data in Fig. 4a (m/cr^ vs. ln(l /cr)) in 
order to discriminate between scalar and fermion behavior. The figure shows a 'window' 
in a where the logarithm of triviality is in the numerator of Eq.(3.9). However, it appears 
that only a small range of a is relevant : for large m (greater than .015) the scaling form 
of Eq.(3.9) probably does not apply, and for small m (less than .002) finite size effects are 
probably distorting the 16^ data. We will gain a better perspective when the 8"^ data is 
plotted. The solid line showing logarithmic behavior is just meant to guide the eye — 
systematic fitting procedures are not appropriate here. 

We also measured the susceptibility at criticality. It follows from eq.(2.4) that, 

^ ^ 21 !i/ A (4-6) 
(T^ ln(l/(T) 

so a plot of x"^/"^^ vs. ln(l/cr) should be linear with a positive slope. Fig.5 shows that 
the data is compatiable with this expectation, but the data suffers from large error bars 
which make precise predictions meaningless here. 

The 8^ data in this case will be subject to much larger and interesting finite size 
effects. The data is shown in Table IV for an estimate of /3c ~ .6066 which follows from 
the equation of state and susceptibilities. Masses ranging from .001 through .020 were 
simulated, and the results are plotted in Fig. 4b in the form suggested by Eq.(3.9). We 
find that the plot ofm/a'^ vs. ln(l./(j) has a negative slope! This should be compared with 
Fig. 4a, where an identical plot was made on the larger lattice (16"^). The two plots differ 
qualitatively. In fact, if we plot a^/m vs. ln(l./a") as suggested by Eq.(3.8), the scalar 
field prediction, then a reasonable plot shown in Fig. 6 follows. This bizarre situation shows 
that finite size effects distort the true behavior of Eq.(3.9), as found on the 16"^ lattice, into 
something compatiable with the scalar field expectation. This should be a severe warning 
for other simulation studies : finite size effects can even obscure the fundamental fermion 
and/or scalar character of the underlying theory. 
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The Average Action and the Specific Heat 

Since accurate measurements of the average action and its associated specific heat are 
relatively easy, let us discuss them briefly. The data is taken from Table I, and the average 
action is plotted vs. (3 in Fig. 7a. We see a clear change in slope near (3 ~ .60. This effect is 
clearer in the specific heat which is plotted in Fig. 8a. As discussed in the theory sections 
above, mean field theory predicts a jump in the specific heat at the critical point. Possible 
curvature and logarithms in such plots are not visible given the error bars in the figures 
and the fact that we could not simulate the model closer to Pc due to finite size effects 
such as vacuum tunnelling. Nonetheless, the (chiral) phase transition is apparent in the 
plots. 

The 8^ data and plots for the average action and its associated specific heat are 
given in the table II and Fig. 7b and 8b. These results are qualitatively similar to the 
16^ results discussed above. They were not useful in subtle studies of logarithms. The 
general agreement between the data from the two lattice sizes is not unexpected since these 
quantities are dominated by ultraviolet fiuctuations. 

6. Discussion and conclusions 

The main conclusions of our study are the following. Logarithmic corrections to 
mean field scaling in the EOS are properly predicted by the large-A?" expansion. Our data 
demonstrate this with a high degree of confidence. In fact, when we tried to impose the 
only remaining possibility, power law scaling in the form of eq.(4.2), not only did the fit 
fail, but it did so in an illuminating fashion - the parameters of the fit flowed toward a 
polynomial generator for the logarithm itself. A pure mean fleld form worked only if a 
restricted range of data were used. 

We have seen that flnite size effects appear in a peculiar fashion and we explained 
the criteria that can detect the finite volume distortions when they are present. These 
criteria are of broader signiflcance since they originate from studies of the simplest of the 
models in the class - NJL with discrete chiral symmetry. Other models with dynamical 
symmetry breaking, like NJL with continuous chiral symmetry or massless fermionic QED, 
are subject to more severe flnite volume effects because of their Goldstone bosons. In that 
context, the most striking feature of our study can be summarized in the comparison of the 
flgures 4a and 4b. They depict identical plots on two different lattice sizes and, because 
of the way the plots are designed, the difference is seen as qualitative. Theoretical input. 
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layed out in ref[20] and summarized in Sec. 3, would have enabled us to discard the 8^ data 
even if we had had no access to the 16^ lattices. While in some instances, the 8^ data were 
not severely distorted (e.g. Fig. 3b), in other instances, they were completely useless even 
as a guide for qualitative behavior. In fact, the logarithms and finite size effects in these 
cases interfere with each other and the data is misleading. Thus, on a small lattice, some 
fits would favor the logarithms in places where they simply can not appear. 

As it turned out, establishing the log's in the NJL model was easier than the corre- 
sponding task in 0^ theory. It was not so crucial to have terribly high statistics as long as 
the lattice size was sufficiently large. Our preliminary results obtained with lower statistics 
did not differ considerably from the final ones. This advantage could be correlated with 
our control over the quality of data and the magnitude of the finite size effects through 
the bound on the exponent 5. 
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Figure captions 

la. The EOS fit Eq.(4.1) to the 16^ data. 

lb. The EOS fit Eq.(4.3) to the 8^ simulation data. 

2a. Same data as in Fig. la fit to the logarithmic EOS Eq.(2.7a). 

2b. 8^ data plotted following Fig. 2a. 

3a. Susceptibility data on the 16^ lattice and power law fits following Eq.(4.4). 
3b. Susceptibility data and power law fits on S'* lattice following Fig. 3a. 
4a. 16^ simulation data at criticality plotted following Eq.(3.9). 

4b. 8^ simulation data at criticality plotted following Eq.(3.9). Note the qualitatively 
different behavior from 4a. 

5. Susceptibility data at criticality plotted as in Eq.(4.6). 

6. 8^ simulation data at criticality plotted following scalar field theory. Finite size effects 
have distorted the data dramatically. 

7a. Average action plotted vs. coupling on the 16^ lattice. 

7b. 8^ analog of Fig.7a. 

8a. Specific heat associated with Fig. 7a. 

8b. 8^ analog of Fig.8a. 
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Table I 

Final compilation 16"^ 



p 


S 


C 




X 


swee 


0.30 


4.8173(1) 


3.8(1) 


1.2123(1) 


0.50(1) 


2 


0.35 


4.5661(3) 


3.6(2) 


1.0084(1) 


0.47(2) 


2 


0.40 


4.3255(1) 


3.44(8) 


0.82828(5) 


0.55(1) 


3 


0.45 


4.0967(4) 


3.438(4) 


0.6613(2) 


0.697(6) 


3 


0.50 


3.8820(4) 


3.11(15) 


0.4974(2) 


0.988(41) 


3 


0.54 


3.7233(3) 


3.01(6) 


0.3592(1) 


1.58(5) 


5 


0.545 


3.7047(2) 


3.01(6) 


0.3408(1) 


1.71(4) 


5 


0.55 


3.6864(2) 


3.01(6) 


0.3220(1) 


1.90(7) 


5 


0.555 


3.6681(1) 


2.84(7) 


0.3025(1) 


2.03(7) 


6 


0.56 


3.6504(1) 


2.80(7) 


0.2826(1) 


2.24(8) 


6 


0.565 


3.6329(2) 


2.79(4) 


0.2619(2) 


2.56(4) 


9 


0.57 


3.6159(1) 


2.81(6) 


0.2403(2) 


3.03(8) 


11 


0.575 


3.5995(2) 


2.79(7) 


0.2179(2) 


3.61(10) 


10 


0.58 


3.5839(2) 


2.73(6) 


0.1945(2) 


4.29(12) 


13 


0.585 


3.5688(1) 


2.59(5) 


0.1693(2) 


5.26(11) 


22 


0.59 


3.5544(2) 


2.55(4) 


0.1415(3) 


7.18(18) 


20 


0.595 


3.5408(2) 


2.66(13) 


0.1091(4) 


12.07(37) 


32 


0.61 


3.5200(1) 


1.002(8) 




24.59(1.76) 


13 


0.615 


3.5194(1) 


0.971(12) 




14.33(96) 


10 


U.OZ 


3.51901(1) 0.968(11) 






Q 
O 


0.625 


3.5186(1) 


0.972(5) 




8.02(37) 


9 


0.63 


3.5184(1) 


0.978(8) 




6.70(23) 


10 


0.635 


3.5181(1) 


0.985(11) 




5.61(18) 


11 


0.64 


3.5178(1) 


0.996(10) 




4.93(17) 


11 


0.645 


3.5176(1) 


1.022(12) 




4.18(17) 


9 


0.65 


3.5173(1) 


1.023(12) 




3.82(13) 


9 
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Table II 

Final compilation S'^ 



p 


S 


C 


a 


X 


swee 


0.54 


3.7385(4) 


2.97(1) 


0.3724(7) 


1.42(5) 


5 


0.55 


3.7036(4) 


2.97(1) 


0.3384(7) 


1.65(5) 


5 


0.56 


3.6699(6) 


2.88(9) 


0.3030(8) 


1.97(8) 


5 


0.57 


3.6367(2) 


2.75(3) 


0.2650(2) 


2.38(3) 


30 


0.575 


3.6206(2) 


2.76(4) 


0.2448(2) 


2.78(4) 


30 


0.58 


3.6052(3) 


2.69(3) 


0.2242(3) 


3.23(5) 


30 


0.585 


3.5890(2) 


2.71(3) 


0.2004(3) 


4.30(12) 


30 


0.59 


3.5739(4) 


2.69(5) 


0.1750(8) 


5.93(29) 


30 


0.595 


3.5584(6) 


2.57(2) 


0.1447(13) 


8.37(9) 


30 


0.62 


3.5234(1) 


1.060(6) 




12.82(19) 


30 


U.ozo 


o.OzzUI i 1 


U.yoD 1 o 1 






oU 


0.63 


3.5208(1) 


0.942(7) 




7.83(19) 


30 


0.635 


3.5200(1) 


0.942(32) 




6.42(11) 


30 


0.64 


3.5193(1) 


0.932(34) 




5.40(7) 


30 


0.645 


3.5188(1) 


0.886(3) 




4.66(7) 


9 


0.65 


3.5184(1) 


0.872(3) 




4.14(10) 


30 


0.655 


3.5180(1) 


0.866(4) 




3.71(10) 


30 
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Table III 

Final Compilation at Criticality [Pc !» 0.6034, 16'') 



m 




X 


sweeps ( 


0.001 


0.0962(3) 


8.87(30) 


22 


0.002 


0.1212(2) 


5.87(4) 


15 


0.003 


0.1389(2) 


4.58(13) 


16 


0.004 


0.1532(2) 


3.68(9) 


16 


0.005 


0.1653(2) 


3.20(8) 


12 


0.010 


0.2093(3) 


2.11(2) 


5 


0.015 


0.2401(1) 


1.68(7) 


5 


0.020 


0.2644(1) 


1.40(2) 


5 
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Table IV 

Final Compilation at Criticality (/5c ~ 0.6055, 8^) 



X 



sweeps (X 



0.1007(9 
0.1360(7 
0.1559(6 
0.1698(5 
0.1808(5 
0.1905(4 
0.1998(4 
0.2076(3 
0.2147(3 
0.2218(3 
0.2278(2 
0.2341(3 
0.2398(2 
0.2457(4 
0.2506(4 
0.2555(3 
0.2594(2 
0.2639(2 
0.2683(2 
0.2726(2 



16.3(8) 
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